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Using data collected at the Pierre Auger Observatory during the past 3.7 years, 
we demonstrated a correlation between the arrival directions of cosmic rays 
with energy above ~ 6 x 10^^ electron volts and the positions of active galac- 
tic nuclei (AGN) lying within ~ 75 megaparsecs. We rejected the hypothesis 
of an isotropic distribution of these cosmic rays with at least a 99% confidence 
level from a prescribed a priori test. The correlation we observed is compati- 
ble with the hypothesis that the highest energy particles originate from nearby 
extragalactic sources whose flux has not been substantially reduced by inter- 
action with the cosmic background radiation. AGN or objects having a similar 
spatial distribution are possible sources. 

Cosmic rays are energetic particles and nuclei from space that strike the Earth's atmosphere. 

Their energies vary from a few 10^ eV to beyond 10^° eV. The flux of cosmic rays at Earth 

decreases very rapidly with energy, from a few particles per square centimeter per second in the 

low-energy region to less than one particle per square kilometer per century above 10^° eV. The 

identification of the sources of ultrahigh-energy cosmic rays (UHECR) with energies ~ 10^° eV 
*The full list of authors and their affiliations appears at the end of this paper. 
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has been a great challenge since they were first observed in 1962 (|7]). Because cosmic rays at 
these energies are not expected to be confined by magnetic fields in the disk of our galaxy, and 
indeed no significant excess from the direction of the Milky Way has been observed, it is likely 
that they originate outside the Galaxy. Until now there has been no experimental confirmation 
of this hypothesis. 

Because of their very low flux, UHECR can only be detected through their interaction 
with the Earth's atmosphere, producing a cascade of billions of particles that excite nitrogen 
molecules in the air along their path and spread over a large area when they reach the ground. 
The Pierre Auger Southern Observatory (|2]), now nearing completion in Argentina, was de- 
signed to simultaneously observe the shower particles at ground level and the associated flu- 
orescence light generated in the atmosphere. A large array of 1600 surface detectors (SDs), 
laid out as an equilateral triangular grid with 1500-m spacing, covers an area of 3000 km^ and 
detects the particles at ground level by means of the Cherenkov radiation they produce in wa- 
ter. At each of four sites on the periphery of the instrumented area, six inward-facing optical 
telescopes observe the sky on clear moonless nights. These devices measure the atmospheric 
fluorescence light produced as an extensive air shower passes through the field of view. The 
two techniques - the SDs and the fluorescence detectors (FDs) - are complementary, and also 
provide cross-checks and redundancy in the measurement of air shower parameters. The SD 
measures the two-dimensional lateral structure of the shower at ground level, whereas the FD 
records the longitudinal profile of the shower during its development through the atmosphere. 
In Figure [TJ we present the layout of the Observatory as of 30 September 2007. 

The Pierre Auger Southern Observatory has been taking data stably since January 2004. 
The large exposure of its ground array, combined with accurate energy and arrival direction 
measurements, calibrated and verified from the hybrid operation with the fluorescence detectors, 
provides an opportunity to explore the spatial correlation between cosmic rays and their sources 
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in the sky. 

If cosmic rays with the highest energies are predominantly protons or nuclei, only sources 
closer than about 200 Mpc from Earth can contribute appreciably to the observed flux above 
60 EeV (1 EeV = 10^*^ eV). Protons or nuclei with energies above 60 EeV interact with the 
cosmic microwave background (|5]|4]|5]), leading to a strong attenuation of their flux from distant 
sources. This attenuation is known as the Greisen, Zatsepin and Kuzmin (GZK) effect, from the 
names of the three physicists that predicted it. If the sources of the most energetic cosmic rays 
are relatively nearby and are not uniformly distributed, then an anisotropic arrival distribution is 
expected, provided the particles have a sufficiently small charge and a sufficiently high energy 
for their directions to be minimally perturbed by intervening magnetic fields. 

Anisotropy of the cosmic rays with the highest energies could manifest as clustering of 
events from individual point sources or through the correlation of arrival directions with a 
collection of astronomical objects. The Akeno Giant Air Shower Array (AGASA) Collabo- 
ration claimed some excess of clustering at small angular scales compared to isotropic expec- 
tations @, but this was not supported by data recorded by the HiRes experiment Q. Analyses 
of data recorded by several air- shower experiments revealed a general correlation with the di- 
rection of the supergalactic plane (|^|9]), where several nearby galaxies cluster, but with limited 
statistical significance. 

AGN have long been considered sites where energetic-particle production might take place 
and where protons and heavier nuclei could be accelerated up to the highest energies yet mea- 



sured {10 11 1. Here, we report the observation of a correlation between the arrival directions 
of the cosmic rays with highest energies measured by the Pierre Auger Observatory and the 
positions of nearby AGN from the 12*^ edition of the catalog of quasars and active nuclei by 
Veron-Cetty and Veron (V-C catalog) ([72]). 
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Data set and method The data set analyzed here consists of the cosmic -ray events recorded 
by the surface array of the Observatory from 1 January 2004 to 31 August 2007. It contains 
81 events with reconstructed energies above 40 EeV and zenith angles smaller than 60°. The 
integrated exposure is 9.0 x 10'^ km^ sr year. 

We only use recorded events if they meet strict criteria with regard to the quality of the 
reconstruction of their energy and direction. The selection of those events is done via a quality 



trigger (75 1 which is only a function of the topology of the footprint of the event on the ground. 
This trigger requires that the detector with the highest signal must be surrounded by five ac- 
tive nearest neighbors, and that the reconstructed shower core be inside an active equilateral 
triangle of detectors. This represents an efficient quality cut while guaranteeing that no crucial 
information is missed for the shower reconstruction. 

The arrival direction of a cosmic ray is a crucial ingredient in our study. The event direction 
is determined by a fit of the arrival times of the shower front at the SD. The precision achieved 
in the arrival direction depends on the clock resolution of each detector and on the fluctua- 
tions in the time of arrival of the first particle (H). The angular resolution is defined as the 
angular aperture around an arrival direction of cosmic rays within which 68% of the showers 
are reconstructed. This resolution has been verified experimentally with events for which two 
independent geometrical reconstructions can be performed. The first test uses hybrid events, 
which are measured simultaneously by the SD and the FD; the second one uses events falling 
in a special region of our array where two surface stations are laid in pairs 11m apart at each 
position. Events that triggered at least six surface stations have energies above 10 EeV and an 



angular resolution better than 1° (15 16). 

The energy of each event is determined in a two-step procedure. The shower size S, at a ref- 
erence distance and zenith angle, is calculated from the signal detected in each surface station 
and then converted to energy with a linear calibration curve based on the fluorescence telescope 
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measurements ([77). The uncertainty resulting from the adjustment of the shower size, the con- 
version to a reference angle, the fluctuation from shower to shower, and the calibration curve 
amounts to about 18%. The absolute energy scale is given by the fluorescence measurements 
and has a systematic uncertainty of 22% ( fT^l ). The largest systematic uncertainty arises primar- 
ily from an incomplete knowledge of the yield of photons from the fluorescence of atmospheric 
nitrogen (14%), the telescope calibration (9.5%) and the reconstruction procedure (10%). Ad- 
ditional uncertainty in the energy scale for the set of high-energy events used in the present 
analysis is due to the relatively low statistics available for calibration in this energy range. 

Events with energy above 3 EeV are recorded with nearly 100% efficiency over the area 
covered by the surface array. The nonuniformity of the exposure in right ascension is below 1%, 
negligible in the context of the present analysis. The dependence of the exposure on declination 
is calculated from the latitude of the detector and the full acceptance for showers up to 60° 
zenith angle. 

A key element of our study is the probability P for a set of events from an isotropic flux 
to contain k or more events at a maximum angular distance from any member of a collection 
of candidate point sources. P is given by the cumulative binomial distribution Y^jLk C^p^{l — 
p)^^^, where the parameter p is the fraction of the sky (weighted by the exposure) defined by 
the regions at angular separation less than ip from the selected sources. 

We analyze the degree of correlation of our data with the directions of AGN referenced in 



the V-C catalog (72). This catalog does not contain all existing AGN and is not an unbiased 
statistical sample of them. This is not an obstacle to demonstrating the existence of anisotropics 
but may affect our ability to identify the cosmic-ray sources unambiguously. The catalog con- 
tains 694 active galaxies with redshifts z < 0.024, corresponding to distances D smaller than 



100 Mpc (19). At larger distances, and around the Galactic plane, the catalog is increasingly 
incomplete. 
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Exploration and confirmation Using data acquired between 1 January 2004 and 26 May 
2006, we scanned for the minimum of P in the three-dimensional parameter space defined by 
maximum angular separations 'ip, maximum redshifts z^ax, and energy thresholds Eth- The 
lower limit for the scan in corresponds to the angular resolution of the surface array. Our 
scan in energy threshold and maximum distance was motivated by the assumption that cosmic 
rays with the highest energies are the ones that are least deflected by intervening magnetic 
fields and that have the smallest probability of arrival from very distant sources due to the GZK 
effect 

We found a minimum of P for the parameters tp = 3.1°, Zmax = 0.018 {Dmax = 75 Mpc), 
and Eth = 56 EeV. For these values, 12 events among 15 correlate with the selected AGN, 
whereas only 3.2 were expected by chance if the flux were isotropic. This observation motivated 
the definition of a test to validate the result with an independent data set, with parameters spec- 



ified a priori, as is required by the Auger source and anisotropy search methodology (20 21 ). 

The Auger search protocol was designed as a sequence of tests to be applied after the ob- 
servation of each new event with energy above 56 EeV. The total probability of incorrectly 
rejecting the isotropy hypothesis along the sequence was set to a maximum of 1%. The param- 
eters for the prescribed test were chosen as those, given above, that led to the minimum of P in 
the exploratory scan. The probability of a chance correlation at the chosen angular scale of a 
single cosmic ray with the selected astronomical objects is p = 0.21 if the flux were isotropic. 
The test was applied to data collected between 27 May 2006 and 31 August 2007, with exactly 
the same reconstruction algorithms, energy calibration, and quality cuts for event selection as in 
the exploratory scan. In these independent data, there are 13 events with energy above 56 EeV, 
of which 8 have arrival directions closer than 3.1° from the positions of AGN less than 75 Mpc 
away, with 2.7 expected on average. The probability that this configuration would occur by 
chance if the flux were isotropic is 1.7 x 10^'^. Following our search protocol and based on the 
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independent data set alone, we reject the hypothesis of isotropy in the distribution of the arrival 
directions of cosmic rays with the highest energies with at least a 99% confidence level. 

Results Having determined that an anisotropy exists, based on the a priori prescription, we 
rescanned the full data set from 1 January 2004 to 31 August 2007 using the method described 
above to substantiate the observed correlation. We used steps of 0.1° in ijj, in the range 1° < 
■0 < 8°, and 0.001 in Zmax, in the range < Zmax < 0.024. We also used a newer version of our 
reconstruction and calibration algorithm that gives slightly different reconstructed directions 
and energies. These small differences, well within our reconstruction uncertainty, modify the 
final event selection, but this has minor consequences on the value of the parameters ip, Zmax, 
and Eth that maximize the correlation signal. We start the scan with the event of highest energy 
and add events one by one in order of decreasing energy, down to Eth = 40 EeV. 

Strong correlation signals occur for energy thresholds around 60 EeV and several combina- 
tions of the other parameters in the range ip < 6°, and Zmax < 0.024 {D^ax < 100 Mpc). The 
absolute minimum value of P occurs for the 27 events with the highest energies (above 57 EeV 
in the new analysis). We generated simulated sets of directions, drawn from an isotropic distri- 
bution in proportion to the relative exposure of the observatory. Performing an identical scan 
on those simulated samples to that applied to the real data, we obtain smaller or equal values of 
P in ~ 10^^ of the simulated direction sets. 

We present (Figure |2]) a sky map in Galactic coordinates of our 27 highest-energy events 
(E > 57 EeV), as determined by our most recent version of the reconstruction code. The 
anisotropy is clearly visible. We note the proximity of several events close to the supergalactic 
plane, and also that two events arrive within 3° of Centaurus A, one of the closest AGN, marked 
in white on the figure. 
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Discussion With the statistics of our present data set, the observed correlation is significant 
for maximum distances to AGN of up to 100 Mpc, for maximum angular separations of up 
to 6°, and for energy thresholds around 60 EeV. Those numbers are to be taken as indicative 
because the minimization of P is not totally exempt from biases. Accidental correlation with 
foreground AGN different from the actual sources may induce bias towards smaller maximum 
source distances while accidental correlation with distant background ones may reduce the op- 
timal maximum angular separation by a few degrees. 

Under the simplifying assumptions of a uniform distribution of sources with equal intrinsic 
luminosity and continuous energy loss in the cosmic microwave background due to the GZK 
effect 90% of the protons arriving at Earth with energy exceeding 60 EeV originate from 
sources closer than 200 Mpc. This (somewhat arbitrarily defined) "GZK horizon" decreases 
rapidly with increasing energy and drops to 90 Mpc for energies exceeding 80 EeV. The relation 
between the horizon distance and the value of Dmax that minimizes P is not a simple one, 
given the possible biases in the method, which has nonuniform sensitivity over the range of 
parameters scanned. Licreasing catalog incompleteness also prevents confidently scanning over 
sources at distances much larger than 100 Mpc. Moreover, the local density and luminosities 
of sources could have significant departures from the uniformity assumed in the GZK horizon 
scale for a given energy threshold. Taking into consideration these caveats, in addition to the 
uncertainty in the reconstructed energies, the range of Dmax and Eth over which we observe 
a significant correlation is compatible with the frequently made assumption that the highest 
energy cosmic rays are protons experiencing predicted GZK energy losses. We note that the 
correlation increases abruptly at the energy threshold of 57 EeV, which coincides with the point 
on the energy spectrum recently reported from the observatory at which the flux is reduced by 
~ 50% with respect to a power law extrapolation of lower-energy observations p7[ ). 

If the regular component of the galactic magnetic field is coherent over scales of 1 kpc with a 
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strength of a few /iG, as indicated by data from studies of pulsars (22), the observed correlation 
over an angular scale of only a few degrees for ~ 60 EeV is indicative that most of the 
primaries are not heavy nuclei. 

These features are compatible with the interpretation that the correlation we observe is evi- 
dence for the GZK effect and the hypothesis that the highest-energy cosmic rays reaching Earth 
are mostly protons from nearby sources. 

The catalog of AGN that we use is increasingly incomplete near the Galactic plane, where 
extinction from dust in the Milky Way reduces the sensitivity of observations. Deflections from 
the Galactic magnetic field are also expected to be significantly larger than average for cosmic 
rays that arrive at equatorial Galactic latitudes, because they traverse a longer distance across 
any regular Galactic magnetic component. These effects are likely to have some impact upon 
the estimate of the strength of the correlation. Six out of the eight events that do not correlate 
with AGN positions within our prescribed parameters and reconstruction code lie less than 12° 
away from the Galactic plane. 

Despite its strength, the correlation that we observe with nearby AGN from the V-C cata- 
log cannot be used alone as a proof that AGN are the sources. Other sources, as long as their 
distribution within the GZK horizon is sufficiently similar to that of the AGN, could lead to 
a significant correlation between the arrival directions of cosmic rays and the AGN positions. 
Such correlations are under investigation in particular for the Infra-Red Astronomical Satellite 
(IRAS) galaxies. The autocorrelation signal of the highest-energy events is also being investi- 



gated. It shows departures from isotropic expectations at angular scales between 5° and 20° (23 1 
and serves as an additional tool to identify the spatial distribution of the sources. 

Conclusion We have demonstrated the anisotropy of the arrival directions of the highest- 
energy cosmic rays and their extragalactic origin. Our observations are consistent with the 
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hypothesis that the rapid decrease of flux measured by the Pierre Auger Observatory above 
60 EeV is due to the GZK effect and that most of the cosmic rays reaching Earth in that energy 
range are protons from nearby astrophysical sources, either AGN or other objects with a similar 
spatial distribution. 

The number of high-energy cosmic -ray events recorded so far by the Pierre Auger Observa- 
tory and analysed in this work corresponds to 1.2 years of operation of the complete southern 
array. The data set that the observatory will gather in just a few more years should offer a better 
chance to unambiguously identify the sources. The pattern of correlations of cosmic-ray events 
with their sources could also assist in determining the properties of the intervening magnetic- 
field structures and in particle physics explorations at the largest energies. Astronomy based on 
cosmic rays with the highest energies opens a new window on the nearby universe. 
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Figure 1 : Layout of the Pierre Auger Southern Observatory. The dots represent the position of 
each of the 1600 SD stations. The 1430 SD stations deployed and activated as of 30 September 
2007 lie in the area shaded blue. The 4 FD sites are labeled in yellow, with green lines indicating 
the field of view of the six telescopes at each site. To give the scale of the Observatory, the 
lengths of the green line correspond to 20 km. 
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Figure 2: Aitoff projection of the celestial sphere in galactic coordinates with circles of radius 
3.1° centered at the arrival directions of the 27 cosmic rays with highest energy detected by 
the Pierre Auger Observatory. The positions of the 472 AGN (318 in the field of view of the 
Observatory) with redshift z < 0.018 (D < 75 Mpc) from the 12*^^ edition of the catalog 
of quasars and active nuclei (12) are indicated by red asterisks. The solid line represents the 
border of the field of view (zenith angles smaller than 60°). Darker color indicates larger relative 
exposure. Each colored band has equal integrated exposure. The dashed line is the supergalactic 
plane. Centaurus A, one of our closest AGN, is marked in white. 
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